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I. INTRODUCTION 



Supergravity was originally developed as an elementary field theory which should avoid 
the ultraviolet divergencies and consequently would represent the long awaited unification 
of gravity with the remaining fundamental interactions of nature. Supersymmetry plays 
an important role in the development of unification models beyond the standard model of 
elementary particles, in the formulation of the conceptual fundamentals of quantum field 
theory and quantum gravity, and more recently in the understanding of important aspects of 
superstring theory (see, for instance, [| for reviews of the conceptual basis of supersymmetry 
and supergravity). Today, supergravity is considered in the first place as an effective field 
theory which should describe the low-mass degrees of freedom of a more fundamental theory, 
probably the still unknown M-theory; however, at the moment the only known candidate 
for such a theory is superstring theory (see, for example, P). 

The classical field equations following from the N — 1 supergravity Lagrangian were de- 
rived in by using the Hamiltonian formalism. There are constraints for each of the gauge 
symmetries contained in the theory: spacetime diffeomorphisms, local Lorentz invariance, 
and supersymmetry. One important result that follows from the analysis of the field equa- 
tions is that one of the constraints relates the torsion tensor and the Rarita-Schwinger field 
so that it can be used to eliminate the torsion tensor from the theory. 

The canonical quantization of supergravity is performed in general by applying Dirac's 
procedure for constrained systems. One uses the 3+1 decomposition of the canonical theory 
to obtain a Hamiltonian in which the symmetry generators of the gauge fields are constrained 
by Lagrange multipliers. Then it is postulated that the wave function is annihilated by all the 
constraints. In the case of N = 1 supergravity there are three constraints: the Hamiltonian 
constraint, the generators of Lorentz rotations and the supersymmetric constraint. It turns 
out jj| that the Hamiltonian constraint is identically satisfied once the supersymmetric 
constraint is fulfilled. Accordingly, only the Lorentz and supersymmetric constraints are the 
central issue. 

The study of N = 1 supergravity models has been limited so far to minisuperspace models 
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12j as direct generalizations of quantum cosmology models. The stan- 



dard approach to quantum cosmology consists in canonically quantizing a minisuperspace 
model which is obtained by imposing certain symmetry conditions on the metrics allowed 
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on the spacelike slices of the universe. This procedure, however, reduces the number of 
degrees of freedom to a finite number and the problem of quantization can be attacked by 
applying the canonical methods of quantum mechanics Q] . The dynamics of the system 
is governed by the Wheeler-DeWitt equation which is a second order differential constraint 
equation following from general covariance, and acts on the wave function of the universe. 
The most general minisuperspace models analyzed in the literature correspond to homoge- 
neous and anisotropic Bianchi cosmological models. Since the corresponding metrics depend 
only on time, the dynamics of the spacelike 3-dimensional slices becomes trivial, unless an 
additional reparametrization is performed. Usually, in the reparametrization one of the scale 
factors of the Bianchi metric is taken as "internal time" so that the Wheeler-DeWitt equation 
generates a wave function of the universe which explicitly depends on the internal time and 
the remaining scale factors. Although this is a quite elegant procedure which in each case 
leads to an explicit wave function of the universe, the main problem regarding the existence 
of classical initial singularity remains unsolved. In fact, in all analyzed minisuperspaces the 
classical singularity remains at the quantum level. Moreover, the original hope that the 
behavior of minisuperspace models would hold at least qualitatively in the full theory seems 
to be not realized. In fact, in Q| it was shown that even in the simple case of a microsuper- 
space (a reduced minisuperspace) contained in the seed minisuperspace the behavior of the 
corresponding wave functions is widely different. The problem of the initial classical singu- 
larity has been attacked alternatively by proposing a wave function for the ground state ly] 
or a tunneling effect Q], among other proposals. In both cases the main idea consists in 
replacing the classical singularity by an ad hoc postulated universe. The initial singularity 
problem remains thus unsolved. The consideration of an additional Rarita-Schwinger field 
in the context of supergravity minisuperspace models (supersymmetric quantum cosmology) 
does not solve the problem, and the classical singularity remains. 

In all the cases the failure to solve the singularity problem can be attributed to the fact 
that, due to the strong symmetry reduction, only a finite number of degrees of freedom can be 
considered. To face this difficulty one needs to analyze genuine field theories with an infinite 
number of degrees of freedom. An option would be to consider milder symmetry reductions 
which leave unaffected a specific set of true local degrees of freedom. These are the so called 
midisuperspace models. Such spacetimes have a long history in general relativity. Indeed, 
any spacetime which allows the existence of two commuting Killing vector fields leads to a 
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real field theory with an infinite number of degrees of freedom. The field equations in this 
case can be shown to be equivalent to the wave equation for a scalar field propagating in a 
fictitious flat 2+1- dimensional spacetime Q]. The local degrees of freedom are contained 
in the scalar field. 



In this work we wi 
mological models 
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1 consider the specific midisuperspace described by Gowdy T 3 cos- 



20j in the context of N = 1 supergravity. We will find an explicit 
expression for the wave function of the universe and will show that the mere consideration 
of a genuine field theory leads to a solution of the singularity problem. In fact, the singular 
behavior of Gowdy models at a certain time of their evolution has been investigated in detail 
at the classical level. It has been shown that the behavior of the metric near the singularity 
corresponds to the so called "asymptotically velocity term dominated" (AVTD) behavior 
(see, for instance j^], for a recent review). We will see that the AVTD singular behavior 
disappears at the level of the corresponding wave function of the universe. 

This paper is organized as follows. In Section |n] we revise the canonical formulation 
of supergravity N — 1 and analyze the Lorentz constraint, following closely notations and 
conventions of j^J . In Section IIHI we present the Gowdy T 3 cosmological models and their 
main properties. Section HVl is devoted to the investigation of the supersymmetric constraint 
and the solutions for the wave function of the universe. Finally, Section IS contains several 
final remarks and the conclusions with indications about different possibilities of generalizing 
the results derived in this work. 

We adopt the following conventions and notations. Indices related to world coordinates 
are denoted by Greek letters. The ones from the middle of the alphabet, i. e. /i, z/, run 
over 0,1,2,3 and the ones from the beginning of the alphabet, i. e. a, /?,..., represent only 
spatial coordinates 1,2,3. Capital Latin indices, i.e. A, B, ... can take the values 0,1,2,3 and 
the small ones run over 1,2,3; both of them refer to to components in an orthonormal local 
frame, for which we use the local metric t]ab = diag(+, — , — , — ). 

II. CANONICAL FORMULATION OF N = 1 SUPERGRAVITY 

The fields of N — 1 supergravity in 4 dimensions are the vierbein e^ A and the Rarita- 
Schwinger gravitino ip^, which is a vector of Majorana spinors. The corresponding La- 
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grangian is given by 

c = \v zi qR - ^ A ^ A7 57^VV , (1) 

where 

D u = d u + -u vAB a AB , (2) 

is the covariant derivative with the Lorentz generators a AB = (1/4)(7 A 7' B — 7 B 7 A ). For the 
■y A matrices we use the real Majorana representation 

n ( a 2 \ , ( ia 3 \ _ / -a 2 \ _ / -ia 1 \ 

7° = I , 7 1 = 3 , 7 2 = I , 7 3 = I, 3 

ya 2 I \ ia 3 J ya 2 J ^0 -ia v J 

in which the anticommutator relation {7" 4 ,7' B } = 2t7 AB is satisfied, and a 1 are the standard 
Pauli matrices. Moreover, 75 = i 7 7 1 7 2 7 3 . In this representation ip = —iip T '-f is known as 
the Majorana condition. The (endomorphic) components of the Ricci rotation coefficients 
can be obtained from Cartan's first structure equation de A = —uj a b A e B , and ojab = 
tOnABdx^. Notice that in general the Ricci rotation coefficients contain a contorsion term. 
However, as we mentioned in the introduction the torsion tensor can be eliminated from the 
theory by using one of the constraint equations. 

Applying the canonical 3+1 decomposition of spacetime j^| the canonical variables can 
be chosen to be the spatial components of the tetrad vectors e a a , their conjugate momenta 
p a a , and the spatial covariant components of the spinor ip a . In this case, the corresponding 
temporal components turn out to be Lagrange multipliers of the Hamiltonian which contains 
only the constraints associated with the three different types of symmetries of the system: 

H = e A U A + ~u AB J AB + S , (4) 

Here 71a contains the usual Hamiltonian and diffeomorphisms constraints, Jab is the Lorentz 
constraint and S denotes the supersymmetric constraint. According to Dirac's canonical 
quantization procedure for constrained systems, the physical states |\&) must be annihilated 
by the corresponding constraint operators, i.e. 

S|tt) = 0, H A \V)=0, Jab\^) = 0. (5) 

From the fact that the supergravity operators satisfy Teitelboim's algebra J^J it follows that 
the condition Ha\^) = is satisfied identically once S\^) = is fulfilled. Consequently, we 
need to consider only the Lorentz and supersymmetric constraints. 



It is convenient to use, instead of the gravitino field, the densitized local components (^o 
is a Lagrange multiplier and so is 4>o ) 



4>a = e e a a ^ a 



as the basic fields commuting with all non-spinor variables, where e 
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det(e a Q ) 



Moreover, if we choose an SO (3) basis it is possible to show that all bosonic terms of the 
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Lorentz constraint cancel each other, yielding |3J 
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and the generator of supersymmetry is given by 



S = e 0al35 l5la D^ s , 

where a factor ordering is usually chosen j^ . 

Let us first analyze the Lorentz condition Jab\^) = that explicitly reads 
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(9) 



Since this constraint does not affect the component a particular solution is to keep 
arbitrary and 

*u = *m = *iv = 0. (10) 

This is the "rest-frame" solution j^J in which only the component \&j remains to be de- 
termined by the supersvmmetric constraint. If none of the conditions is satisfied, the 
Lorentz constraint implies that 



(11) 
(12) 

(13) 



Notice that in the representation we are using here the Lorentz generators J a ^ are 4x4 
matrices. Therefore, the components ^n, ^in, and \I/jy of the wave function must be 
considered as 4 x 1 matrices. 
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There are two alternative ways to solve the system of algebraic equations given in (jllj) - 
(fHIj) . The first one consists in taking the components ^n, ^m, and \1/ jy proportional to 
each other. The proportionality factors can be absorbed by redefining each component of 
the wave function so that this case can be written as 

^ii = ^iii = ^iv ■ (14) 
Then from Eas.([TT |) -([T3 |) it follows that 

>Jl2 = —>Jl3 — <^23 ) (15) 

a condition which implies a trivialization of the Lorentz constraint. Then from the expression 
(JTj) we can determine the components of the gravitino field. In the trivial case (fTo'j) we obtain 



^3 • 



(16) 



The second possibility is to solve explicitly the system of equations (jTTj) - ([T3j) . The rela- 
tion between the components of the Lorentz generators can be solved by representing each 
component as a product of 7-matrices, under the condition that the corresponding Lorentz 
algebra is preserved. Alternatively we can use the standard generators of the ordinary 
rotation group 0(3) as given in |2J] 

/o oo\ /ooo \ 
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Solving explicitly the Lorentz constraint PH ) - (fT3*j) in this representation it is easy to see 
that for each of the vectors ^n, ^111, ^ iv only one component is different from zero. That 
is to say the wave function of the universe can be represented as 

/ \ / a \ 
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^111 
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\d J 



(18) 



where E is a function to be determined, a is an arbitrary 4-vector and b , Co, and do 
are arbitrary constants. For the sake of simplicity later on we will consider the special case 



a o = { fl o, 0, 0, 0}. So we have shown that the wave function (|18|) explicitly solves the Lorentz 
constraint. The remaining function E and the arbitrary constants ao, bo, cq, and do have to 
be chosen so that the supersymmetric constraint is satisfied. 

This ends the analysis of the Lorentz constraint. Notice that in the "rest-frame" solution 
(|1U|). the wave function of the universe is a scalar with only one independent component. For 
the further investigation of this solution as well as of the non-trivial solution (|14j). we need 
to analyze the supersymmetric constraint (jHJ), in which the bosonic part plays a prominent 
role. In the next Section we will present the specific gravitational field which completely 
determines the bosonic part of the supersymmetric constraint. 



III. GOWDY T 3 COSMOLOGICAL MODELS 



Gowdy cosmological models are inhomogeneous time-dependent solutions of Einstein's 
vacuum equations with compact Cauchy spatial hypersurfaces whose topology can be either 
T 3 or S l x S 2 

Other particular topologies are contained in these two as special 
cases. Here we will focus on T 3 models for which the line element can be written as 

ds 2 = e- x ' 2+T ' 2 {e- 2T dT 2 - d X 2 ) - e~ T [e p (da + Qd5) 2 + e~ p d6 2 ] , (19) 

where P, Q, and A depend on the non-ignorable coordinates r and \. The spatial hypersur- 
faces (r =const) are compact if we require that < \i °i $ < 27r. The expression in square 
brackets depicts the metric on the T 2 subspace which is generated by the commuting Killing 
vectors d a and d$. The coordinate \ labels the several tori. 

When the Killing vectors are hypersurface orthogonal, the general line element (|19|) be- 
comes diagonal with Q = and the corresponding cosmological models are called polarized. 
In this last case, the subspace T 2 corresponds to the spatial surfaces of a 2 + 1 fictitious flat 
spacetime in which a scalar field, represented by the metric function P, propagates j3, El- 
The local degrees of freedom contained in the scalar field are true gravitational degrees of 
freedom which cannot be eliminated by a choice of gauge. We are thus facing a genuine field 
theory which is a special case of a midisuperspace model. Notice that the infinite number 
of degrees of freedom contained in this midisuperspace model can be associated with the 
inhomogeneous character of the spacetime. If we would neglect the inhomogeneities present 
in the model, we would obtain a minisuperspace model with a finite number of degrees of 
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freedom, probably related to a Bianchi cosmological model. The general unpolarized case 
(Q 7^ 0) also corresponds to a midisuperspace model; however, its interpretation in terms 
of a dynamical scalar field in a 2+1 spacetime can not be realized. In this work we will 
concentrate on the polarized case (Q = 0) where the field equations can be integrated in 
general. The general unpolarized case (Q ^ 0) will be also considered in quite general terms 
at the level of the wave function of the universe, although no exact solution with Q ^ will 
be analyzed due to the difficulty of the classical field equations. 

The vacuum field equations for the general line element (|19|) can be written as a set of 
two second order differential equations for P and Q 

P TT - e~ 2T P xx - e 2P (Q 2 T - e- 2T Q 2 x ) = , (20) 
Q TT - e~ 2r Q xx + 2(P T Q T - e~ 2T P x Q x ) = , (21) 

and two first order differential equations for A 

A r = P r 2 + e- 2 -P 2 + e 2P (Q 2 T + e~ 2 ^) , (22) 
A x = 2{P X P T + e 2P Q x Q T ) . (23) 

Notice that the function A can be calculated by quadratures once P and Q are known. The 
system of differential equations (J2(J|) and (J21j) for P and Q is highly nonlinear, has been 
investigated in detail by using numerical methods and has been analyzed analytically 



only recently in 
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where several special solutions have been derived. In the special 
polarized case, using the method of separation of variables it is possible to find the general 
solution as 

oo 

Q = 0, P(r, X ) = Y} An cos ( n *) + B » sin [CnMne~ T ) + D n N (n e - T )] , (24) 

n=0 

where A n , B n , C n and D n are arbitrary constants and Jo, No are Bessel functions. The 
integration of the function A is quite cumbersome for the general solution. We will present 
in Section HVl a particular solution characterized by a finite number of terms of the sum (|24j) . 

The behavior of Gowdy cosmological models near the singularity is an important property 
that has been intensively used to study the geometric behavior of the initial Big-Bang 
singularity of our Universe. In the case of T 3 models it can be shown that the singularity is 
approached in the limit r — ► oo. It has been proved that all polarized Gowdy models belong 
to the class of "asymptotically velocity term dominated" (AVTD) solutions and it has been 



conjectured that the general (unpolarized) models are also AVTD 13311 . This conjecture has 
been reinforced through the analogy with other midisuperspaces |3J]. The AVTD behavior 
states that near the singularity each point in space is characterized by a different spatially 
homogeneous cosmology j^]. It implies that at the singularity all spatial derivatives of the 
field equations can be neglected and only the temporal behavior is relevant. This leads to 
a "truncated" set of differential equations which in the case of T 3 models can be obtained 
from Eqs. pUj) - (j2HJ) by neglecting all the derivatives with respect to the spatial coordinate 
X- It is easy to see that the general solution to this "truncated" system is given by j^ll Q 

Pavtd = ln[a(e~ CT + b 2 e CT )] , Qavtd = , ^ , , 2 . + d , Xavtd = A + c 2 r , (25) 

where a, b, c, d, and Ao can be considered as arbitrary real constants. The singularity 
situated at r — > oo is characterized by a blow up of the curvature which is determined by 
the behavior of the AVTD solution (1251) . 

For the calculation of the supersymmetric constraint we need explicitly the connection in 
an orthonormal basis. The structure of the line element (|19|) suggests the following choice 
for the vierbein 

e = e -(A+3r)/4 rfr ^ g l = g -(A-r)/4 dx ? g 2 = e (P-r)/*( d(T + Q dS ) , £ 3 = e ~{P^ dS , ( 26 ) 

which satisfies the orthonormality condition g^e li e B u = r] AB with e A = e^dx^. In this 
tetrad the non-vanishing components of the connection lo^ab are 

^ooi = — ^A x , u; 23 = — —e + Q T , 

w ifii = -^ r ( 1 ~ > ^123 = \ e? Qx » 

^202 = X -e^~^ (1 - P T ) , W 2 3o = -\^Qr , 

1 3_i A , P _ 1 3_, A , 3n„ . , 

^2i2 = ~c^~ X , ^213 = 2^ * 2 Qx > (27) 



W 302 



^303 



\e^ +P [{l ~ e- T P T )Q ~ Qr] 



1 



312 ~~ 2 v xV ' 

^313 = -\e-^-i{P x - e * p Q x Q). 

where the hat refers to indices associated to the local orthonormal basis. 
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IV. PHYSICAL STATES 



In this section we analyze the remaining supersymmetric constraint. The densitized local 
components of the gravitino field (jUJ) depend on the components of the local Gowdy tetrad 
(J2SJ). Noting that in this case e = = det(e a a ) = exp[(A + 3r)/4], we obtain from Eq. (JHJ) 
that 

^ = e -i(A-H-)^ , ^ = e -K A+5 ^0 2 , ^ = e"^ A+5T+2P )(0 3 + e P Q0 2 ) . (28) 

With these values for the gravitino field and the connection components (|27[). it is now 
straightforward to determine the explicit form of the supersymmetric constraint JSJ) which 
after lengthly calculations can be written as 

S = e-fa+^Si + e^ x+2r ^S 2 , (29) 



S 1 = i 7 ° 7 1 



~ \\ + \Px) + 7 3 03 U - \\ x - -P x J + ^ 7 3 2 e P Q x 



(30) 



5 2 = z 7 °{[e~ P (7 2 7 3 03 - 7SV1) + Q(7S 3 0i + 7 2 T 3 02)]9 CT - (7 V02 + l 1 ^^} ■ (31) 

The component Si contains all the terms which include dependence on the non-ignorable 
spatial coordinate \- The second components contains derivatives with respect to the spatial 
coordinates a and 5. Since the classical spacetime metric depends only on the spatial 
coordinate \i one could expect a similar dependence for the wave function |\&). In such a 
case, the action of S2 on the wave function would vanish and we would need to consider 
only the first term S\. However, we will analyze here the most general case allowed by the 
supersymmetric constraint S\^) = 0. 

According to the discussion of Section [HI there are two different types of solutions to 
the Lorentz constraint: the rest-frame (jlOj) . and the trivial solution (|14j). The next step 
is to solve the supersymmetric constraint for each of these special solutions. However, as 



mentioned above the rest-frame solution leads to scalar wave functions. This imp 
the wave function contains only bosonic degrees of freedom and, according to 



ies that 
the 



corresponding states cannot be physical. In fact, physical states must contain fermionic 
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degrees of freedom and have infinite number of modes. In the case of the trivial solution 
(|14[). its non-physical character could have been expected from the fact that the trivialization 
of the Lorentz constraint as given in Eq.(|15|) does not fulfill Teitelboim's algebra ^j. Thus, 
we are left only the non-trivial solution (|18j) . On the one hand, this special non-trivial choice 
of the gravitino field guarantees that fermionic degrees of freedom will enter the final form 
of the wave function and, on the other hand, it satisfies identically the Lorentz constraint in 
the sense that it leads to an expression for the wave function which involves the fermionic 
variable in a manifestly Lorentz invariant combination. 

From the general expression for the supersymmetric operator (jSHJ)? (JHDI), and ({SI} , we 
obtain 



Si 



(e" p r 4 + QT 5 )d a - T% 



4 A x 4 P ^) r2 + 2" i 



e"Q x r 3 , 



(32) 
(33) 



where, in order to solve them, we use the following realization for the product of the 7 
matrices with 

,3 1 



-Z7 



in terms of the following T l matrices 



(34) 



r 1 
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-1 
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-1 
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1 








(35) 



/ 11 \ 
-10 0-1 
-10 0-1 
\ 11 / 



/ -1 1 \ 
10 10 
10 1 
\1 -10/ 



(36) 



The physical states are given as the non-trivial solutions of the equation S\^>) = 0. 



221 showed that there are no 



The investigation of this equation in minisuperspace models 
non-trivial solutions and, therefore, sometimes it is believed that supergravity iV = 1 is 
an uninteresting theory with no physical states. A less radical conclusion would be that 



minisuperspace models are unphysical due to the fact that the strong symmetry reduction, 



12 



which leads to a system with finite number of degrees of freedom, does not allow the existence 
of non-trivial physical solutions to the supersymmetric constraint. In the present work we are 
dealing with a midisuperspace with an infinite number of degrees of freedom and therefore 
the last explanation does not hold. Indeed, in this case the non existence of genuine physical 
states would indicate a very serious difficulty for N = 1 supergravity. We will show that 
this is not the case. 

From the above considerations it is clear that a family of physical states can be obtained 
from wave functions which simultaneously satisfy the constraints <Si|\l/) = and c^l 1 ^) = 0. 
As mentioned above, if in accordance to the functional dependence of the classical metric 
we suppose that the wave function depends only on the spatial coordinate \i the constraint 
$2^) — is identically satisfied and we only need to solve the set of differential equations 
following from «Si|W) = 0. If we limit ourselves to wave functions which are independent 
of the coordinates a and 5, we guarantee that "anomalies" do not appear. Indeed, the 
classical symmetries associated with the Killing vectors (d^) 11 and (ds)^ does not hold at the 
quantum level if we find wave functions which depend on these coordinates. This would be 
an indication of the existence of anomalies at the quantum level. We will show in the next 
subsections that it is not necessary to assume independency of a and 5. We will see that 
the differential equations following from the constraint S\fy) = can be solved by applying 
the method of separation of variables and that the resulting compatibility conditions for the 
wave function of the universe eliminate the possibility of existence of anomalies. 

A. The polarized case 

As mentioned in Section UTTl the polarized case of Gowdy models corresponds to the limit 
Q = of the line element ()19|). and the differential equation for the function P allows the 
general solution given in Eq. ()24|) . The supersymmetric constraint S\^) = leads to the 
following set of first order partial differential equations 




(37) 



d x (V n - vl/,) -\{\ x - P x )^n + \{\ x + P x )*i 
_ e -|(A-2P-3r) [e-PQ^j + q, IV ) + ds ^ !l + vfr m )] 



0. 



(38) 
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_ e -i(A-2P-3r) 



[e- p ^(*/ + + d s {V n + *jv)] = 



(39) 



-<9 x (^ m + ^ 7 y) + i(A x - P x )^ /y + i(A x + P x )* m 



0. 



(40) 



First, we consider the special case where the components of the wave function are inde- 
pendent of the spatial coordinates a and 5. This means that the constraint ^I^P) = is 
identically satisfied and the last set of equations reduce to a system of ordinary differential 
equations. To find solutions to the resulting system it is natural to use the exponential func- 
tions as an ansatz for each of the components of the wave function. It is then straightforward 
to show that the general solution is given by 



I*) 



^ii 
^iii 
V ^iv ) 



b 
("0 
\d / 



(41) 



where ao, bo, Cq, and do are arbitrary constants satisfying the relationships 



og + 6g = 



cg + dg = 



(42) 



for which we choose the solution a — an d c = =R(io- This wave function of the universe 
represents a physical state for the N = 1 supergravity Gowdy model. It is interesting 
to notice that in this particular case we were able to completely integrate the system of 
differential equations and all the components of \^/) are explicitly given in terms of the 
"classical" functions P and A. Moreover, the metric function P enters the wave function of 
the universe in a way very similar to a phase which does not affect the physical significance 
of the solution. Indeed, if we define the "absolute value" of the wave function of the universe 
as |*)|*)* = + ^ n ^* n + ^hi^ui + ^jy^/y, we obtain from Eq.flU]) 



|tf )|#>* = (a a* + Mo + cqCq + do<Qe*< 



(43) 



Notice that this "absolute value" can not be associated with a probability density (see 
the discussion in Section M). The behavior of the wave function of the universe is thus 
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entirely dictated by the metric function A. Let us recall that the function P obeys a linear 
differential equation and also determines the behavior of A through Eqs. (j22|) and (|2*Hj) which 
are nonlinear. We see that the nonlinear sector of the classical field equations enters the 
final form of the wave function of the universe, whereas the linear sector appears as a phase 
only 

With the wave function of the universe we can analyze the problem of the cosmological 
singularity. Recall that the classical spacetime is characterized by a singularity at r — > oo 
where the metric functions behave according to the AVTD solution (|25|). In particular, the 
metric function Xavtd diverges linearly as r — > oo. This singular behavior is illustrated in 
Fig. 01 Let us consider a special solution contained in (|24|) with D n = and C n — 1. We 
choose this special case for the sake of simplicity and because it reproduces the structure of 
the general solution. Consider also the case with n — 1, 2, i.e. 



Notice that the first term n = leads to a constant value of P which can be absorbed in 
the metric through a coordinate transformation and leads to the Minkowski metric. For this 
reason we do not consider this term in the series (j2H). Using the expression (I44|) the field 
equations (|2^|) and (f2l!j) can be integrated and yield 



The wave function of the universe for this special case is depicted in Fig. El where we 
can see that it is regular for all values of \ and r. As the classical singularity (r — > oo) 
is approached, the wave function of the universe remains constant and finite. This result 
shows that the classical cosmological singularity of this particular Gowdy model has been 
removed after its canonical quantization in the context of N = 1 supergravity. The presence 
of the fermionic field is thus sufficient to solve the singularity problem at the quantum 
level in this particular midisuperspace model. To obtain this result we have considered in 
(I44J) only the first two nontrivial terms of an infinite series. If we consider further terms, 




(44) 




15 




FIG. 1: The wave function of the universe is formally defined for the AVTD solution as \^>)avtd ~ 
exp(AAVTD/2) = exp[(Ao + c 2 r)/2], according to Eg. 1)25(1 . The graphic shows the singular behavior 
of \^)avtd, for r — > oo, with c = 1, and Ao = 1. 

the expression for the corresponding function A become rather cumbersome and difficult to 
handle. Nevertheless, since all higher terms contain only Bessel functions, the mathematical 
structure of the function A resembles that of Eq. ()45|) so that for further terms we can expect 
again wave functions of the universe which are free of singularities. To confirm in an invariant 
way that no singularities appear in the wave function of the universe we have analyzed the 
curvature spatial scalar R = e^R a a , which characterizes the structure of the constant time 
hypersurface where the wave function of the universe is denned. In the present case we find 
that R = — (1/2) exp[(A — r)/2]P 2 , and a numerical analysis of this scalar shows that it is 
regular everywhere. 

Let us now consider the general case in which the components of the wave function depend 
also on the spatial coordinates a and 5. Using the method of separation of variables, it is 
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FIG. 2: Behavior of the wave function of the universe for the special solution (|44|l. and (|45|). The 
AVTD singular behavior for r — > oo, has disappeared and it is replaced by a well-behaved function 
with a constant value. 

easy to show that the general solution to the system (j3*Tj) - (j4Tjj) must have the form 



I*} 



VlII 

V Viv ) 



D ma+n5+g(x) 



b 
c 
\d J 



(46) 



where m and n are constants which arise from the separation of variables, and g(x) is a 
function to be determined. Furthermore, ao, b Q , c , and do are arbitrary constants. Inserting 
this ansatz for the wave function into Eqs. (j37|) and (JHHj) . we obtain one single differential 
equation 



1 



i(A-2P-3r) 



if the following relationships are satisfied 

CoOo - ao) + d (a + b ) = , 

Similarly, from Eqs. ()39|) and (J4"U|) we obtain 
1 . i 



-p Oo + ^0 

e m- 



a + c 



&o — ao ~ «o 

d (a - 6 ) + c (a + &o) = 



(A-2P-3t) 



_p bo + Co ao — Co 

e m n — 

Cq + cto Cq + ao 



(47) 



(48) 



(49) 
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when the conditions (|4"H)l are fulfilled. Notice that Eqs.ffS]) imply that eft + b% = and 
Cq + c?q = so that we can choose 

a = ±ib , c = =Rdo , (50) 

as the general solution for the system (|48[). Comparing term by term the two equations for 
the function g x given above we find two new conditions on the constants bo and do, i.e. 

(a + d )(c + d ) = (b + c )(b - a ) , (a + c )(c + d ) = (a - c )(a - &o) • (51) 

Notice that these conditions follow from the terms containing the constants "m" and "n" 
in Eqs. fl47[) and (|49jl so that they are valid only if m ^ and n ^ 0. It is now easy to see 
that Eqs. (|51jl does not allow any solutions compatible with (|5()jl. Consequently, the only 
compatible solution is for m = = n, a result which implies that an explicit dependence on 
the spatial coordinates a and 5 is not allowed and, therefore, the existence of anomalies in 
the wave function of the universe is completely excluded. 

Under the above considerations, the resulting differential equation for g(x) can easily 
be integrated and yields g(x) = (A =F iP)/4. The imaginary part can again be considered 
as a phase and, therefore, the behavior of the wave function of the universe is completely 
dictated by the behavior of the metric function A only. The explicit final expression for 
coincides with the one of the former case given in Eq. (j41j) . 



B. The unpolarized case 

In this Section we will analyze the general unpolarized T 3 Gowdy model (Q ^ 0). No 
general solution to the classical field equations is known in this case because of their highly 
nonlinear character. In fact, only a few exact highly nontrivial solutions are known in 
the H_ H fl Q Q. Let an exantpte be gi ve n for the function, P, Q, and A 
satisfying the field equations (|2()jl - (J2BJ). This is all the information we need to know in 
order to investigate the supersymmetric constraint S\^) = which yields the following set 
of equations 

d x ^j + *„) - \{\ x - P x )9j - \{X X + P x )*„ - l -e p Q x ^ m 

+e -|(A-2P-3r) [ e ~P da ^ H + y m ) _ Qd^jj - * rv ) + d s {*n - = , (52) 
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d x (y 



11 



e -|(A-2P-3r) 



9j) - -(\ x - P x )* n + -(A x + PJVj + -e p Q^ IV 

[e~ p d a (^! + * 7y ) - + ^ui) + + *///)] = , (53) 



_ e -i(A-2P-3r) [ e -^ (j (^ / + Wjy ) _ g^(* /7 + * /V ) + + = , (54) 

-0 x (*m + tf JV ) + ^(A x - P x )^ /y + ^(A x + P x )tt m + -e p Q x V n 
+e -i(A- 2 p-3r) [e^^ + + g 5ff (^ 7 _ _ ds ^ I _ $ m )] = o . (55) 

It is clear that an exponential ansatz with separation of variables similar to the one used 
in the previous case in Eq. (j4fi)) will lead to a system of ordinary differential equations and 
a set of algebraic equations. The analysis of the resulting equations is similar to the one 
performed in the last subsection for the polarized case. In a similar manner, it is possible 
to show that no anomalies are allowed in the wave function of the universe which therefore 
turns out to depend on the spatial coordinate x only. The resulting wave function of the 
universe can be expressed as 



Mx) 



V ®iv J 

where the set of constants must satisfy the conditions 



b 

Cq 

\d J 



(56) 



c (fy) - ao) + d ( a o + b o) = , d (a - b ) + c (a + b ) = 



(57) 



which allow the solution clq = ±ib , cq = =pi<io- For the sake of simplicity we consider the 
case in which b and do are real constants. Then, the function h(x) can be put in the form 



Kx) 



e P Q x dx 



(58) 



1, d f P i f d 

-X-—JeQ x d X = F - i \ < P + - 

This function together with the conditions (|57|) represent an explicit solution of the 
super symmetric constraint and show that in the unpolarized case there exist nontrivial 
physical states. The imaginary part of h(x) can again be interpreted as a phase which does 
not affect the behavior of the wave function of the universe (J23). The behavior of the real 
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part of h(x) is governed by the behavior of the metric function A and an integral which 
involves the functions P and Q. Hence, we need the explicit form of the metric in order to 
analyze the wave function of the universe. As mentioned before, it is very difficult to solve 
the system of differential equations which follow from Einstein's vacuum field equations. To 
analyze the behavior of the wave function near the singularity it would be necessary to apply 
numerical methods for handling the metric functions. This is a task which is outside the 
scope of the present work. 

V. FINAL REMARKS AND CONCLUSIONS 

Our analysis of the supersymmetric constraint presented in the last sections relies on 
a very specific foliation of spacetime. In fact, the general form of the constraint S = 
e 0al3S 'y5 r y a Di3ips fixes at the very beginning the time coordinate x° so that it does not appear 
as a dynamical variable in the further analysis. Here we have chosen x° = r which is 
the natural time coordinate in Gowdy cosmological models. The canonical quantization 
formalism forces us to "freeze" this time during the entire analysis and it reappears explicitly 
only in the solutions for the wave function of the universe, where we interpret it as a label 
which associates a particular value of the wave function to a different spatial slice. It is in 
this sense that we can say that we know explicitly the wave function of the universe at each 
moment of time. And it is in this context that we were able to investigate the behavior of 
the wave function near the cosmological singularity. 

Consequently, the problem of the "frozen" time in our analysis cannot be solved in the 
context of the canonical quantization formalism. The common concern about the use of 
this formalism is whether the final result of the quantization (in our case, the wave function 
of the universe) depends on the choice of a particular foliation. To clarify this question 
in the present case let us consider a different quite general foliation. Consider a new time 
coordinate t defined by 

dt = e-^ x+3T) dr . (59) 
Then, from the general line element (J19j) we obtain 

ds 2 = dt 2 -e 2A d X 2 -e F (e p da 2 + e~ p d5 2 ) , (60) 
where A = A(t,x) and F = F(t, x)- For the sake of simplicity we are considering here 
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the polarized case only (Q = 0). In this particular parametrization the lapse function 
becomes a constant. We apply now the canonical quantization procedure of supergravity 
with a foliation determined by x° = t =const. Clearly, this is a quite radical change of 
foliation when compared with the original one (r =const). To calculate the supersymmetric 
constraint we proceed as in Section HVI for the non-trivial solution (|18|). Then we obtain 

S = e A+2F 5! + e l ^ A+3F+p ^S 2 , (61) 

with 

* = (^x + Ax + ^x + ^^+^ + A.-^ + ^r 2 , (62) 

s 2 = e - p r 4 <9 CT - r 5 d 5 . (63) 

If we compare these expressions with the supersymmetric constraint in the original foliation, 
given in Eqs. (j32j) and (}3*3*|) . we see that the general structure does not change. The integration 
of the differential equations following from applying (|fi2*|) and (jHSJ) to the wave function 
leads to a solution with a structure similar to that of Eq. (f4~Tf . This shows that the general 
behavior of the wave function of the universe is qualitatively invariant with respect to a 
transformation of the time coordinate as given in Eq. (j59|) . This coordinate transformation 
of time is quite general since it involves all the nonignorable coordinates present in Gowdy 
models and implies a constant lapse function. Clearly, the qualitative independence of the 
wave function of the universe holds also in the case of less general coordinate transformations. 

It is necessary to mention that for the expression of the Lorentz constraint in terms of 
the gravitino field ((7J) there is only one possible choice, namely 

0i ~ 7 1 , 02 ~ 7 2 , 03 ~ 7 3 • (64) 

On the other hand, for the solution of the supersymmetric constraint we use the choice 
which is more convenient for the analysis of the resulting differential equations. Nevertheless, 
one could use the representation ()64j) since the final result for the wave function of the 
universe does not depend on the particular used representation. 

It is important to emphasize that the physical interpretation of the wave function of 
the universe |\&) presents certain difficulties. A genuine wave function must be related to 
observable quantities and this implies that must yield a probability density. However, 
this is not true in this case, in particular because the wave function of the universe is 
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not normalizable. Moreover, if we require that yields a probability density for the 3- 
geometry which must have a sp ecific value at a given time, this would imply a violation of 
the Hamiltonian constraint [15[ . These difficulties in the interpretation of the wave function 
of the universe are the price one has to pay for applying the canonical quantization procedure 
which involves the isolation of a specific "time" parameter against which the evolution of 
the system can be defined. An alternative procedure like the Dirac quantization based on 
functional integrals, which does not require to single out the time variable, could lead to a 
quantum system with less interpretation difficulties 

In this work we have investigated the Gowdy T 3 cosmological models in the context 
of N = 1 supergravity. The quantum constraints resulting from the canonical quantization 
formalism were explicitly analyzed, and for the resulting set of differential equations we were 
able to find general solutions. In this way, we found the wave function of the universe for the 
polarized and unpolarized special cases. This represents a proof of the existence of physical 
states in the (N = 1) supersymmetric midisuperspace corresponding to Gowdy cosmologies. 
This result contrasts drastically with analogous investigations in minisuperspace (Bianchi) 
models where no physical states exist, a result that sometimes is assumed as a sufficient proof 
to dismiss supergravity N = 1. We have adopted a less radical position in this work and 
dismiss as unphysical only the minisuperspace models. The existence of physical states in 
midisuperspace models confirms this conclusion and indicates that supergravity N = 1 is a 
valuable theory which should be investigated further. In this context we have also obtained 
an interesting result showing that, in the Gowdy T 3 midisuperspace model analyzed in 
this work, the wave function of the universe which represents nontrivial physical states is 
completely free of anomalies. 

Gowdy cosmologies are characterized by the existence of a curvature singularity at tem- 
poral infinity, and it is known that the metric functions near the singularity evolve according 
to the AVTD behavior. We have shown that after the canonical quantization the correspond- 
ing wave function of the universe is free of singularities. This represents a solution to the 
singularity problem and is one of the main results of this work. The mere presence of the 
Rarita-Schwinger field and the consideration of a genuine midisuperspace model is sufficient 
to eliminate the classical singularity. This result points to a further interesting and expected 
property of supergravity N = 1 in the sense that it is able to properly handle the conceptual 
limits of classical general relativity. 
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In this work we focused on the special case of T 3 cosmologies. The generalization of 
our results to include the case of S 1 x S 2 Gowdy models seems to be straightforward. In 
particular, we believe that the unified parametrization introduced in (3?}], which contains 
both types of topologies, could be useful to explore the supersymmetric Gowdy model in 
quite general terms. 
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